Josephson current through a quantum dot coupled to a molecular magnet 
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Josephson currents are carried by sharp Andreev states within the superconducting energy gap. 
We theoreticaUy study the electronic transport ol a magnetically tunable nanoscale junction consist- 
ing of a quantum dot connected to two superconducting leads and coupled to the spin of a molecular 
magnet. The exchange interaction between the molecular magnet and the quantum dot modifies 
the Andreev states due to a spin-dependent renormalization of the quantum dot's energy level and 
the induction of spin-flips. A magnetic field applied to the central region of the quantum dot and 
the molecular magnet further tunes the Josephson current and starts a precession of the molecu- 
lar magnet's spin. We use a non-equilibrium Green's function approach to evaluate the transport 
properties of the junction. Our calculations reveal that the energy level of the dot, the magnetic 
field and the exchange interaction between the molecular magnet and the electrons occupying the 
energy level of the quantum dot can trigger transitions from a to a tt state of the Josephson 
junction. The redistribution of the occupied states induced by the magnetic field strongly modifies 
the current-phase relation. The critical current exhibits a sharp increase as a function of either the 
energy level of the dot, the magnetic field or the exchange interaction. 

PACS numbers: 74.45. +c, 74.50.+r, 73.23.-b, 75.50.Xx 
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I. Introduction 

Molecular spintronics combines the two fields molecu- 
lar electronics and spintronics.'^' In spintronics, the elec- 
tron spin is used as the degree of freedom in which in- 
formation is encoded.'^ Molecular electronics investigates 
the electrical and thermal properties of molecules and 
aims to build devices composed of single molecules or en- 
sembles of molecules. In molecular spintronics, the spin 
of molecules is used to manipulate the spin and charge 
transport. In particular, molecular magnets are interest- 
ing as basic building blocks for electronic devices^^^HI and 
for quantum computing.'^ ® These molecules have a per- 
manent magnetisation due the their anisotropy barrier 
as well as long coherence timea^ that facilitate further 
quantum-mechanical phenomena such as int erferenc e^'^ 
and quantum tunnelling of the magnetisation.EESEi] 

Experimentally, the transport properties of different 
kinds of junctions containing magnetic molecules have 
been extensively studied in three-terminal devices .'^^HEl 
Current measurements through molecular magnets al- 
lows to identify the magnetic states and directly observe 
the magnetic anisotropy and the orientation of the easy 
axis.-'^- These magnetic states of the molecule have been 
proposed to enable quantum computing.'SI An alterna- 
tive way to probe the properties of molecular magnets 
is to deposit the molecules on carbon nanotubea^ or 
graphene layers .'^° The presence of the magnetic molecule 
modifies the transport through the junction and oppo- 
sitely, the tunnelling electrons modify the magnetisa- 
tion of the molecule and can reverse the magnetisation.'^GI 
Other experiments used superconducting electrodes with 
the advantage that the heat losses in these devices dis- 
appear. The proximity-induced superconductivity and 
the accompanied Andreev reflection s mod ify the trans- 
port properties through the molecule.'^iEll jjj a,n Andreev 
reflection process, incoming electron-(hole-)like quasipar- 



ticles with energies lying within the superconducting gap 
are retroreflected as hole- (electron-) like quasiparticles at 
a normal-superconductor interface. In a junction con- 
sisting of two superconducting leads coupled over a non- 
superconducting region, the Andreev reflected electron- 
and hole-like quasiparticles at the left and right interfaces 
form Andreev levels which carry the Josephson current. 
Josephson junctions offer the possibility for applications 
in superconducting electronics, as well as quantum in- 
formation and computing. ^"^ ^'^ The Andreev levels of a 
Josep hson junction can in principle be used as a two level 
qubit.'^^'^^'ln Ref. [17] j an Andreev level qubit with spin 
orbit coupling was discussed as building block to perform 
quantum computations using the spin degree of freedom 
in order to manipulate the Andreev states. The manipu- 
lation of the Josephson current by adding quasiparticles 
to the Andreev states and the resulting suppression of 
the Josephson current were measured in Ref. [58 . An 
alternative way to manipulate the Andreev states of a 
quantum point contact in the presence of a magnetic 
scatterer was studied in Ref. [29]. The presence of a 
molecular magnet in a Josephson junction offers a fur- 
ther possibility to tune the current and manipulate the 
Andreev states. These states are experimentally accessi- 
ble and have been observed in Ref. [30]. By comparison 
of the energy of the Andreev states with the experimen- 
tal data it is also possible to extract information about 
the parameters affecting the current in the constriction. 

Besides the spectroscopy of the Andreev states, the 
measurement of the Josephson current reveals detailed 
information about the internal structure of the junc- 
tion. Direct measurement of the current-phase relation 
in superconducting atomic contacts have been performed 
in Ref. [31]. The shape of the current-phase relation 
strongly depends on the details of the contact between 
the electrodes and is important for applications in super- 



conducting devices. The state of the system can change 
from a to a TT state in which the current changes sign. 
This transition to the tt state was proposed in Ref. [52] as 
a result of tunnelhng through magnetic impurities. Ex- 
perimentally, this transition was measured in a junction 
consisting of two superconducting leads coupled over a 
quantum dot in the Couloumb blockade regime which was 
occupied with a odd number of electrons.'^-'' The singlet 
wave function of the Cooper pairs experiences a phase 
shift of TT due to coherent cotunnelling processes, leading 
to the reversal of the supercurrent. The to tt transition 
has been studied in many other kinds of Josep hson j unc- 
tions, such as ferromagnetic heterostructures lM135| q^^^^ 
theoretically analysed in magnetic junctions associated 
with a molecular magnet in Refs. 



In this paper, we study the Josephson current through 
a quantum dot which is coupled to a molecular magnet. 
We focus on the limit of low temperature (T — )• 0) and 
on the regime of negligible Coulomb interaction. The 
magnetic moment of the molecular magnet is assumed 
to be large enough to allow for a classical treatment 
of the molecular magnet's magnetisation. Furthermore, 
we assume that the molecule has an isotropic magneti- 
sation and may, for instance, b e a fu llerene molecule 
doped with a magnetic impurity.'^IESl ^pj^g gpjj-^ Qf ^j^g 

magnetic molecule then interacts with the electrons oc- 
cupying the quantum dot via the exchange interaction. 
In Refs. [35] and [SZ], the Josephson current through an 
isotropic magnetic molecule was studied in the Kondo 
regime and in the regime of negligible Coulomb interac- 
tion, respectively. In comparison to the work in Ref. ^36J , 
in this paper a magnetic field is applied to the central 
region consisting of the quantum dot and the molecu- 
lar magnet, whose magnetisation then precesses with the 
Larmor frequency, lol- The magnetic field and a gate 
voltage applied to the quantum dot introduce additional 
parameters for manipulating the Andreev states and tun- 
ing the current through the junction. We mainly focus on 
the dependence of the Josephson current on the energy 
level of the dot and the magnetic field. The manipulation 
of the Andreev states by the molecular magnet is studied 
and the parameter range enabling the junction to be in 
a or TT state can be determined by the critical current. 



The outline of the paper is as follows. In Sec. [TT| 
we introduce the model Hamiltonian of the junction and 
describe the terms including the effect of the molecular 
magnet. The following Sec. 



Ill is concerned with the 



approach used to determine the transport properties of 
the junction. The results, starting with the density of 
states of the quantum dot and the Andreev states, are 
presented in Sec. 



IV 



This section also includes a dis- 
cussion about the current-phase relation and the critical 
current. We conclude with a summary of the results in 
Sec. El 



II. Model 

The junction, which is depicted in Fig. fll (a), consists 
of two superconducting leads coupled to a quantum dot 
in the presence of a molecular magnet. A magnetic field 
applied along the z axis induces a Zeeman splitting of 
the quantum dot's energy level as well as a precession of 
the spin of the molecular magnet. The Hamiltonian of 
the junction is written as 

Hit)= J2 {Hc.+HTc.) + HD + Hsit) + HsD{t), a) 

a=L,R 

where the index a corresponds to the left [L) or right 
(i?) side of the junction. The left and right supercon- 
ducting leads are described by the BCS Hamiltonian 

He = Efe„^6„4„,TCfcc<T + A„4^.^c^_j^^^ + Al,c_fe^^Cfe„t 
with a — (t, i) = ±1 and the dispersion 
^k — ^ k"^ /{2m) — ji. The order parameter is given by 
and for symmetry reasons we assume 



A, 



lA, 



o^Va 



that ipH.L = i'/'/S. In general, the temperature and 
magnetic field dependence of the order parameter must 
be taken into account in a self-consistent way. Since 
we restrict our discussion to the limit T — > 0, we 
neglect the temperature dependence and for simplicity 
we assume that the applied magnetic field does not 
fundamentally affect the superconducting leads. Quasi- 
particles in the leads with momentum k and spin a 
are created and annihilated by the operators cj.^ and 
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FIG. 1. Junction with two superconducting leads and a quan- 
tum dot coupled to the spin S of a molecular magnet. In the 
laboratory frame (panel (a)), the magnetic field forming an 
angle of j? with the spin, starts a precession of the spin and 
splits the energy level of the quantum dot into Eoi^ and Eo^. 
In the frame of the rotating spin (panel (b)), the Fermi ener- 
gies of the spin-up (spin-down) electrons and the quasiparticle 
states in the leads are shifted by —{+)liJl/'2. The transforma- 
tion cancels the Zeeman splitting of the quantum-dot energy 
level induced by the magnetic field. 



Cka- The Hamiltonian describing the tunnelhng between 
the dot and the superconducting leads is written as 

Htcc = J2k^a '^l^d„:ka.Ck^<y + cl^<T^L,k^^<y The hopping 
yd„,ka describes the couphng between the dot and the 
leads and is assumed to be independent of energy. 
The operators acting on the dot correspond to d\ and 
da- The Hamiltonian of the dot Hd and the Zeeman 
energy are written as Hd — J2cri-^o + criiBBz)d\d„ 
and Hs{t) = — /x^f(i)B, with the magnetic field 
B = (0,0, -B^) and the magnetic moment of the molecu- 
lar magnet fii^f{t). The magnetic moment is related to 
the spin of the molecular magnet via fJ'Mi^) = ~l^i^) 
with the gyromagnetic ratio 7 = qm^/ {2m). The Lande 
factor of the molecular magnet is denoted by gu and 
has in principle to be determined by comparison with 
experiments. In the following, we assume that the Lande 
factor of the molecular magnet equals to that of free 
electrons and gM = 2.^^ Due to the magnetic field, the 
magnetisation of the molecular magnet precesses with 
the Larmor frequency. We assume that the motion of 
the spin is un damped which can be achieved by dc and 
rf fields.HZHll xhe equation of motion is then given by 
dS/dt — — 7SXB. The solution of this equation is S(i) = 
S (cos(cL'Li)sin(i?)e2; + sin(a;Li)sin(^)ej^ + cos(^)ez), 
with the magnitude of the spin |S| ^ S and the Lar- 
mor frequency ujl = jBz- The exchange interaction 
between the molecular magnet and the quantum dot is 
described by Hsoit) = ^J2acr' '^sdl{S{t)(T)„a'da' with 
the coupling Vs between the spin and the quantum dot 
and the Pauli matrices cr — {(Jx,(Jy,(7z). The term can 
be transformed into Hsoit) = J2a'^''^sCos{-d)d'lda + 
v,sin{^)e-"^^*d\di + v,sin{^)e"^^^dld-^ with v^ = SVJ2. 
The first term induces a spin-dependent shift of the 
energy levels of the dot while the second and third terms 
account for the spin flip of the electrons occupying the 
dot. Since the magnetic field enters in the Hamiltonian of 
the dot and the exchange interaction, we can rewrite the 
Hamiltonian of the dot in terms of the Larmor frequency 
of the molecular magnet as Hd = J^a {-^0 + '^^) d^d^. 

III. Approach 

The transport properties of the system are described 
by a nonequilibrium Green's function approach.**^ ^^ In 
order to simplify the evaluation of the Green's func- 
tions, we perform a unitary transformation to the ro- 
tating frame of the molecular magnet's spin, since in 
this frame the Hamiltonian is time independent. The 
state vector transforms according to \^) — W\^) and 
the Hamiltonian in the rotating frame can be written as 
H — U^ HU + i{dtU^)U with the unitary transformation 
operator 



C/(t)=exp 



2 



Olka 



In the rotating frame, the Hamilton operator (nj) is 
given by H=J2JH^+Htc.)+Hd+Hs+Hsd- The 



transformation results in a spin-dependent shift of 
the quasiparticles' energies in the leads and the 
quantum dot. The Hamiltonian of the leads reduces 



Ho.= Ek^.i^k^--'^^)claCk^a+^C 



k^t -k^i 



to 

A^c^kai'^kat- ^^^ ^o ^he spin-dependent energy 
shift, the Zeeman energy in the Hamiltonian of 
the dot is cancelled by the transformation such 
that the Hamiltonian of the dot is given by 
ffjj = J2cr ^odlrda- The spin is fixed in the rotat- 
ing frame and the exchange Hamiltonian is written as 
HsD= J2cr crVsCos{d)dld„+Vsam{i3)d\di+VsSm{d)dld^. 
The remaining terms of the Hamiltonian M are not 
affected by the transformation (l2|. Fig. fl](b) depicts the 
system in the frame of the rotating spin. The splitting of 
the spin-up and spin-down energy levels of the quantum 
dot appears because of the exchange interaction and the 
energy shift of the spin-up (spin-down) quasiparticles is 
given by -f (— )wsCOs(??). 



In order to evaluate the Green's functions of the sys- 
tem, we artificially divide the structure into three subsys- 
tems, which are the left lead (L), the right lead (R) and 
the quantum dot (D).ISSl The spin dependence and the 
superconducting state are taken into account by writing 
the Green's functions in Nambu-spin space. Since the 
system is out of equilibrium, we additionally write the 
Green's functions in Keldysh space. In Keldysh-Nambu- 
spin space, the Green's functions have the structure 



Gi3p'{t,t') 






(3) 



The symbols " and " denote a matrix in Keldysh-Nambu- 
spin and Nambu-spin space, respectively. The labels i?, 
A and K indicate the retarded, advanced and Keldysh 
element of the Green's function Gi3p'{t,t'). The indices 
/? and /3' refer to operators in one of the three subsys- 
tems, e.g. the retarded Green's function G^^ is given 

by G^a{t,t')= - t6{t - t'){{Mt),i^Ut')}), where we 
have introduced the operators of the quantum dot and 
the leads in Nambu-spin space as ^d = [d^ d^ dl dV)'^ 

and ^pa — {ca^ c^i c_^^ '^-ai)'^ with the index a re- 
ferring to the momentum k^ in the left or right lead. 
The elements of the matrix in (Is]) are related by G^ = 
(1/2) (G^ -G^ + G^) with the lesser Green's function 

defined by G<^{t,t') = i(4(t')^i5(i)). 



Taking into account all k states of the leads, we 
define the matrices Gppr and V^^/ as {Gi3j3i)ij=Gp.^i3i. 

and {y i}fi')ij—Vp.,pi, with the diagonal matrices 

Va,d^<iia.g{yd^^a,Vd^,a,-Vd^-a,-Vd^-a) in Nambu- 
spin space. The indices i and j indicate all k states in 
the leads. The matrices (Is]) of all subsystems are then 
combined in an enlarged Hilbert space into one matrix 
defined by G. The full and the unperturbed Green's 



function are written as 

/Gll Old Glr\ ^ /Gol 

G = I Gdl Gdd Gdr. I I Go = I Gqd 
\GflL Grd GuiiJ \ GoRj 

The coupling of the quantum dot to the leads and to the 
molecular magnet is given by 

Vii3 \ 

^ DL ^ DD ^ DR I I 
^RD / 

where Vdd = Ws (<5'zC0s(??) + (T2;Sin(z9)), is treated as a 
perturbation. The elements containing the coupling be- 
tween the quantum dot and the leads are diagonal in 
Keldysh space and are in Nambu-space given by Va^d- 

Fourier transforming the Dyson equation to energy 
space, we obtain 



G = Go + GoVG 



(4) 



and can calculate the retarded and advanced Green's 
functions of the dot as 






yR/A 



C^nn : (5) 



where the self-energy is^ defined as S = VoLgLLVLD + 
VoRgRRVRD = f^L + ^R- These self-energies are ob- 
tained by summation over all quasiparticle states in the 
left and right leads, respectively. This summation can 
be replaced by the integral ^j,^ — ?► (2^)3 ! d^^^ 



-^ 



Nq J d£,a J ^^° with the normal density of states at the 
Fermi energy, A^o- The integration over the quasiparticle 
energies in the superconduc tor le ads to the quasiclassical 
Green's function denoted atP^'^ 

MA ^ -^ ( E""^ A^iay \ 

^" y|A|2 _ (iji?M)2 \tayAl ~E"/^) ' 

together with the normalisation condition iga ) — 

— TT^l. The retarded and advanced self-energies are ap- 
proximated by Ea' — Na\V\'^ga and the energy 
j^R/A jg defined as i?^/^ = E ± i-q with 77 -j' 0. In 
the rotating frame, the self-energies are given by (setting 



E'^+Oi 



^ 



^|A|2-(BH+^)2 







V|Ap-^(£;«-LD)2 V|A|2-(B«-cD)2 











A+, „ „ -{e"+Cj) 



\V|AP-(B«+cD)2 











V|A|2-(£;«+cD)2/ 



(6) 

with the tunnelling rates defined as F^ — nNoVf^^. 
The electron and hole part of the spin-up and spin- 
down unperturbed Green's function of the dot in 



Nambu-spin space is given by {Gqjj)^\^^—E^—Ec, and 

(G'oi))33/44=^(-^^+-^o)- The unperturbed Green's func- 
tion and the self-energy enable us to evaluate the full 
Green's function of the dot ^ . 

The average charge current operator in the Nambu- 
spin space from the left lead to the quantum dot is ob- 
tained by using the Heisenberg equation of motion 



Jl 



H,Nl 



Ht,Nl 



(7) 



with the Hamilton operator H in the rotating frame, the 



2 ^k 



/,t 



and the tun- 

+i^lylkyk^- 



number operator N[^ — 

nelling operator Hra = 5 J2k^ ^jlj^d,k^i^d 
The Josephson current ([7]) can then be written in terms 
of the lesser Green's function as^ 
dE 



Jl- 



2n / 27r 



Tr (<7o {G<lYld-VdlG<o)) . (8) 



By using the Dyson equation Q in the enlarged Hilbert 
space, we calculate the elements G^^ and G^^^ with the 
help of the relation G< = (1/2) (G^ - G" + R^). The 
results are^Sl 



Gld 



G< ^' "^ 



'Dd'^DLSlL + Gno'^DLgLL 

Sll^ldG^^ + Sll^ldGI^^ 

The Josephson current then simplifies to 
dE 



Jl- 



2-K 



Re 



Trao Gg^E 



R v< 
DD^L 



^DD^L 



(9) 
(10) 

(11) 






obtained from the 



G%ot<G'^ij with t< 



The Green's function 

Keldysh equation by G^^ 

Ti^F—FT,^. Due to the transformation to the rotating 

frame, the effective Fermi energy in Nambu-spin space 

are shifted and F is given by 

/f{E+u) 

f{E~Cj) 

f{E~Cb) 





F = 




(12) 



with the Fermi function f{E) = 1/{1 + cxp{E/kBT)). 

IV. Results 

The charge transport properties of the junction can 
now be investigated. In principle, the current is given 
by two contributions. The first is the current carried 
by Andreev states whose energies lie within the super- 
conducting gap. The second contribution is carried by 
continuum states outside the superconducting gap. 



A. Density of states of the quantum dot 

Figure |2] shows the spin-resolved density of states of 
the quantum dot at ip = 7r/2, _Bo = and a symmet- 
ric coupling to the leads with F = F/, = F^ = 0.1 A. 
In panel (a), we consider a static magnetisation of the 



molecular magnet in the z direction (wl = 0, i^ = 0). 
The exchange interaction between the molecular magnet 
and the quantum dot lifts the spin degeneracy and shifts 
the energy level of the quantum dot in the rotating frame 
by ±Ws for spin- up and spin-down electrons, respectively. 
Therefore, the spin-up energy level, which is broadened 
by r, is closer to the upper superconducting gap edge in- 
creasing the density of the spin-up continuum states for 
_B > A whereas the density of the spin-down continuum 
states is decreased. The exchange interaction also lifts 
the spin degeneracy of the Andreev states and shifts the 
states El and En (Em and Eiy ) of the spin-up (spin- 
down) quasiparticles to higher (lower) energies. For the 
parameters in panel (a), the exchange coupling pushes 
both spin-up (spin-down) states above (below) i? = 0. 

In panel (b), a magnetic field is applied in addition 
to the exchange coupling of the molecular magnet and 
the quantum dot. The magnetic field shifts the contin- 
uum states by — (-f)aji/2 for spin-up (spin-down) elec- 
trons due to the transformation in the rotating frame and 
also slightly pushes the Andreev states towards |£'| — > 
compared to panel (a). If the magnetisation points in 
an arbitrary direction (t? ^ 0) , the electrons can undergo 
spin flips into sidebands separated by the energy Ulol ac- 
cording to the Hamiltonian (II]) . In the frame of the rotat- 
ing spin, the transformation compensates the exchange 
of energy and the electrons are scattered into states at 
the same energy. In panel (c), ?? = 7r/4 and the Andreev 
states as well as the continuum states can be occupied 
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FIG. 2. (Color online) Spin-resolved density of states of the 
quantum dot. In (a) and (b), Vs — 0.3 A and i} — 0. In (b) a 
magnetic field of oji, — A/2 is applied, shifting the continuum 
states of the spin-up (spin-down) electrons by ±ij;i/2. In 
(c) and (d), i? 7^ and an electron on the quantum dot can 
change its spin direction such that the states of the spin-up 
and spin-down electrons are mixed. In (d), the magnetisation 
is precessing in the plane {iS = 7r/2). The other parameters in 
panel (a)-(d) are tp = n/2, T = O.IA, Eg = and r; = IQ-^A. 



with spin-up and spin-down electrons respectively. In 
panel (d), '& = 7r/2, the magnetisation is precessing in 
the xy plane and the density of the scattered states in- 
creases. 

The complete parameter dependence of the Andreev 
states in the central region are obtained by the poles of 
the Green's function in Eq. ^ which are given by 

A+A- -B = 0. (13) 

Using the notation cj = lul/2, the elements are 



A±^ Uet VsCos^)^/A^ - {E±Cb)^ + 2r{E ± Cj) 
- 4r2A2cos2((^/2) - E^ (A2 -{E± Cu)^) 



and 



B = w2sin2,?^(A2 -{E- c:;)2)(A2 - {E + tZ))2) 

+ VA^ -{E- (I')2^A2 -{E + cD)2 {-vlsin^i) 
+2 {E^ + Ei - vjcos^'d)) + 

4t({E + Q){E + i;^cos??) ^A2 - {E - Cb)^ 

+{E -Co){E~ w,cos?9)^A2- {E + Oj)A 



The pole equation (13) reduces to that of Ref. [56], if 



the quantum dot does not interact with the molecular 
magnet and no magnetic field is applied {vs = w = 0). In 
this case A+ is equal to A- and S = 0. If no magnetic 
field is applied, the i^-dependence of the Andreev states 
vanishes since no spin-quantisation axis is preferred. In 
this limit, the equation of the Andreev states agrees with 
the result in Ref. [36]. If 1? = 7r/2, the spin precesses in 
the xy plane and the equation of the Andreev states is 
symmetric with respect to lol -^ —^l- Additionally, the 
Andreev states are symmetric under the transformation 

Eq — )• —Eq. 

In the limit £■ ^ A, the equation of the Andreev states 
can be explicitly solved with the result 



E 



1 



iji 



{2TAcos{ip/2)y 



A-^-Cj-^ 



VsCos{'d)~ 



Tlj 



yA2 



+ w2sin^(^) (14) 



and 



^ 1 (2rAcos(^/2)r 
En I, IV - - ( ±V ^2^-2 +^0 



UsCos(z9)- 



Tuj 



VA2-w2 



+ v^,sm\d) (15) 




FIG. 3. (Color online) Phase dependence of the Andreev 
states. In panels (a)-(c), F — A/10 and Eq = 0. In (a), 
the exchange coupling Vs splits the Andreev state into spin- 
up and spin-down states. The efTect of a finite i9 is shown in 
panel (c). Panel (d) shows the Andreev states for F = A/2, 
§ = 7r/4 and a position of the energy level on the dot of _Bo = 
and Eo = A/2. The dashed lines indicate the effective Fermi 
energies of the spin-up and spin-down quasiparticles at zero 
temperature. 



with V = y|Ap"^^^/ (y^\A\^~^'^ + 2rj . In general, 

however, the Andreev states must be calculated numer- 
ically from Eq. ( 13 ) and are shown in Fig. pj In panel 



(a), the spin degeneracy is lifted due to the exchange 
interaction between the quantum dot and the molecular 
magnet. The spin- up (spin-down) electrons are shifted to 
higher (lower) energies. Since the Andreev states lie well 
inside the energy gap, we can use Eqs. (14 1 and (151 to 



find expressions for the splitting of the Xndreev states. 
The displacement of the Andreev states in panel (a) due 
to the exchange interaction is given by ±Vs/{A + 2r). 
The exchange interaction shifts the Andreev states across 
the Fermi energy, which from equation ( 12 ) is located at 
_E = at r = 0. Therefore, the current is expected to 
be strongly modified in panel (a) if the coupling v^ is in- 
creased. An applied magnetic field counteracts the shift 
of the Andreev states induced by the exchange interac- 
tion, Vs- The combined shift of the Andreev states due 
to Vs and ujl is given by ±{l/v){vs — Tu}/{y/A'^—uj'^)). 
The effective Fermi energies in panel (b) at T = and 
uiL = A/2 are located at the energies —{+)ujl/'2 for spin- 
up (spin-down) electrons. In this case, both spin-up An- 
dreev states are shifted above the effective Fermi energy 
Ep^, whereas both spin-down Andreev states are below 
the effective Fermi energy Ep^. The effect of ■& on the 
Andreev states is shown in panel (c). Now, all Andreev 
states below Epi (Ep^) are occupied with spin-down (up) 
electrons similarly to the situation of the density of states 
in Fig. [2] (c) and (d). In (d), F = A/2 and the shift of 
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FIG. 4. Current-phase relations for F = A/10 and _Bo = 
(panel (a)-(c)). In (a) lol ~ and j? = and Va varies from 
to 0.3A. The values correspond to the Andreev states of 
Fig. [3] (a), where the current is strongly suppressed due to 
the shift of the Andreev states across the Fermi energy. In 
panel (b), Vs — 0.3A, •& = and ujl varies from to A. In 
(c), Vs = 0.3 A, ljl ~ A/2 and & varies from to n/2. In (d), 
Vs = 0.3A, LOL = A/2, ^ = 7r/4, F = A/2 and Eo is varied. 
The temperature is set to ksT — 10~*A and -q — 10~'*A. 



the Andreev states due to F is larger than the splitting 
due to Vs such that at ip = the states Ejj and Ejy are 
below the effective Fermi energy Ep^ and the states Ej 
and Ejj J are above the effective Fermi energy Ep^. 



B. Current-phase relations 

The Andreev states discussed in the last section carry 
the Josephson current /n(<^)- The contribution to the 
current of each Andreev state En{ip) is proportional to 
the derivative of the Andreev state with respect to the 
phase multiplied with the Fermi function at the energy 
of the Andreev state,!^ 



n dip 



(16) 



The Josephson current is then given by the summation 
over all Andreev states n. 

In the following, we discuss the current-phase rela- 
tion obtained by Eq. (11) in the low temperature limit 
ksT — 10~^A. This expression contains the contribu- 
tion to the current from Andreev as well as the contin- 
uum states. Fig. H (a) shows the current-phase relation 
of a static spin in the z direction {■& = 0) and zero mag- 
netic field {ojl = 0) and four different values of Vs ■ The 
Andreev states corresponding to the values of Ws =0 




/(eA/^i) 0.2 



FIG. 5. Current-phase relation for Vs = 0.1 A, i9 — 7r/4, ljl ~ A/2 (a). By changing Eq, 
■K state. In (b), So = Q,Vs = 0.3A and i9 = QAn. The magnetic field changes from uil ~ 
r = Tl = Tfl = O.IA, ksT = IQ-^A and rj = 10~''A. 



the state changes from a to a 
— A to cjl = A. In both panels 



and Vs — 0.3A are depicted in Fig. |3] (a). At Vg — 0, 
the spin-degenerate Andreev state below the Fermi en- 
ergy Ep carries a positive current for (p < t: according 
to relation (16 1. If the exchange coupling is increased to 
Vs — 0.3A, the contributions from the current of both 
current-carrying Andreev states Ejjj and Ejv in Fig. [3] 
(a) cancel whereas the Andreev states Ej and Ejj are 
completely unoccupied above the Fermi energy. In this 
case, the current is strongly suppressed since in total the 
Andreev states do not contribute to the current and the 
current is carried by the continuum states giving rise to 
a TT state of the junction. Between Vs = and the com- 
plete shift of the Em state below the Fermi energy, the 
current-phase relation sharply decreases at phases where 
the Em level intersects the Fermi energy. These phases 
are determined from Eq. (13) to ip — 2arccos(tis/2r). 
At phases larger than 2arccos(tis/2r), a negative current 
appears due to the continuum states since the Andreev- 
state contributions to the current cancel. 

Figure [4] (b) shows the current-phase relation for the 
same parameters as in panel (a) at Vg = 0.3A and 
i} — 0, but the magnetic field is increased from uj^ = 
to ujl = A. The Andreev states corresponding to the 
current-phase relation at uj^ = A/2 and wi, — A are de- 
picted in Fig. [3](b). For the parameters chosen in|4](b), 
the spin-up (spin-down) Andreev states do not cross the 
corresponding effective Fermi energy Ep^ {Epi) and no 
sharp change of the current is observed. The Andreev 
states El and En [Fig. [3](b)], which carry spin-up elec- 
trons, are for all values of a;^ in Fig. |4](b) above the effec- 
tive Fermi energy Ep^ and therefore do not contribute to 
the current. The Andreev states Em and Eiv carrying 
spin-down electrons are both below the effective Fermi 
energy Ep\^ and do not contribute either due to current 
cancelation of Em and Eiv state. Therefore, only the 
continuum states give rise to a current such that the junc- 
tion is in the tt state. 

In Fig. H (c), the current-phase relation is shown for 



the same parameters as in panel (b) at lj^ = A/2 but 
the angle d increases from to Tr/2. Since the electrons 
on the quantum dot can undergo spin flips at finite i?, 
the Andreev states in the rotating frame are degenerate 
and we have to take into account four spin-degenerate 
Andreev states following from the density of states in 
Fig. [2] (c) and (d) . The Andreev states corresponding to 
the current-phase relation in Fig. [4] (c) at ^ = 7r/4 and 
I? = 7r/2 are shown in Fig. [31(c). At i? = 0, the current is 
the same as in panel (b) at w^ = A/2 where Ei and En 
are empty while Em and Eiv are occupied with spin- 
down electrons and therefore the Andreev states do not 
contribute to the current. If d is increased, spin-down 
(spin-up) electrons are scattered into Ei and En {Em 
and Eiv) states. The particles scattered into the En and 
Eiv states are below the corresponding effective Fermi 
energy Ep^ and Ep^, respectively. These two states give 
a positive contribution and therefore a finite d increases 
the current. The sharp step for phases close to tt of the 
current-phase relation appears because an Andreev state 
crosses the Fermi energy, similarly as in panel (a). For 
r = A/2, the current-carrying Andreev states cross the 
Fermi energy twice thus leading to the two steps in the 
current-phase relation as shown in panel (d) . A finite Eq 
opens a gap of the Andreev states dX Lp — t:. 

So far we have discussed the current-phase relations 
of the Andreev states shown in Fig. [3J We now con- 
sider how the to TT transition and the reverse process 
are driven by Eq or uji^. The current-phase relation as a 
function of Eq is shown in Fig. [s] (a). The state changes 
from a tt junction at £^0 = to a junction at Eq = A. 
The current increases stepwise as a function of p, but the 
position of this stepwise increase changes as a function of 
Eq. This behaviour appears twice, since the spin degen- 
eracy of the Andreev states is lifted due to the exchange 
interaction. Since the current is symmetric under the 
transformation Eq — >■ —£^01 the same transition is driven 
by decreasing Eq from Eq = Q to Eq = —/S.. In panel 



(b) of Fig. [5J the current-phase relation is evaluated as 
a function of lj^ ■ If the magnetic field is increased from 
uji = —A to cji = 0, the junction is driven from a tt 
to and back to a tt state. Since the Andreev states 
are symmetric under the transformation lo^ — >■ — wl if 
1} — O.Stt (Sec. IV A I, the Josephson current shows a 



similar behaviour if ?? approaches O.Stt. 



C. Critical Current 

As we discussed in the previous section, different pa- 
rameters can drive the junction from a to a tt state 
or vice versa. In order to further investigate the trans- 
port properties, we consider the critical current, which is 
experimentally more easily accessible than the current- 
phase relation. 

Figure [6] shows the critical current as a function of 
the energy level at lol = A/2 and different couplings 
Vs in panel (a). For small Eq, the state of the junction 
undergoes a tt to transition by increasing Vg and finally 
goes back to the tt state. If the junction is in a tt state 
at Eq, a transition is driven to the state by increasing 
£^0 for all values of the exchange coupling shown in panel 
(a). The critical current at Vs = 0.1 A in the range of 
Eq — to Eq = A/2 corresponds to the maximal current 
of each current phase relations shown in Fig. ^ (a). The 
sharp increase of the critical current occurs due to the 
shift of Andreev states across the Fermi energy within a 
small parameter range of Eq. In panel (b), Vg — 0.3A and 
the critical current is depicted for different values of lul ■ 
Again, by increasing Eq, the junction undergoes a tt to 
transition at ujl = and oj^ = ±1.2A. In addition, the 
junction also exhibits multiple transitions for Eq < A/2 
as ujl is increased. 

Contour plots of the critical current are depicted in 
Fig. [Tj In all panels, we assume a symmetric coupling 
between the quantum dot and the leads and calculate the 
critical current in the low temperature limit. In (a) and 
(b) , the critical current is plotted as a function of Eq and 
Vs a,t d = 7r/4 and the magnetic field is set to oj^ = in 
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FIG. 6. Critical current as a function of Eq. The parameters 
in (a) corresponds to Fig. ^ (c) for different exchange cou- 
plings at i9 = 7r/4 and ujl ~ A/2. In (b), Vs — 0.3A and 
'& = 7r/4. In both panels F = A/10 and ksT = IQ-^A. 
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FIG, 7, (Color online) Critical current Ic in units of eA/h as 
a function of Eo and t^s at 1? = 7r/4 (panel (a) and (b)). In 
(a) cji = and in (b) ujl ~ A/2, In (c) and (d), the critical 
current is shown as a function of Eq and ljl at Vs = 0,3 A, 
The angle i9 varies from t? = 7r/4 in (c) to i9 = 0,47r in (d). 
The other parameters are F = Fl = Tr — 0,1A and ksT = 
10-*A, 



(a) and A/2 in (b). In panel (a), the junction exhibits 
a to TT transition due to an increase of the exchange 
coupling, Vs- In this case, the contribution to the current 
of the Andreev states does not depend on 1? since the 
magnetic field is zero. The Andreev states intersect the 
Fermi function at Vs = ±^4F2cos((pmax/2)^ -I- Eq, with 
the phase <pmax corresponding to the phase at the critical 
current. At these values of the exchange coupling, the 
junction exhibits the transition, livg > 2F|cos((^,„a2:/2)|, 
the junctions also exhibits tt to transitions as a result 
of increasing Eq. In panel (b), the magnetic field is set 
to UJL = A/2 and the Andreev states are degenerate due 
to the spin-flip scattering. The parameters in this panel 
correspond to the critical current shown in Fig, ^ (a). 
Transitions are driven by either increasing Eq or Vg . Due 
to the redistribution of the occupied states at finite Wi , 
the critical current exhibits an additional step which is 
due to a shift of the Andreev states across the Fermi 
energy. In panel (c) and (d) , the critical current is shown 
as a function of Eq and ujl at Vs = 0.3A and i? changes 
from -d ~ 7r/4 to ?9 = 0.47r. Due to the spin-flip term, the 
critical current shows a mirror structure around ujl — 
with different intensities for ujl — ^ — w/,. As discussed in 
Sec. |IVA[ the Andreev states are symmetric with respect 
to UJL -^ — wl at I? = 0,57r, Therefore, the critical current 
shows a symmetric behaviour as a function of ujl since 
-d is close to O.Btt. For small Eq, the junctions exhibits 
multiple transitions as a function of the magnetic field. 
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V. Conclusions 

We have studied the electronic transport in the low 
temperature limit of a Josephson junction consisting of 
two superconducting leads and a quantum dot connected 
to a molecular magnet. The exchange interaction lifts the 
spin degeneracy of the energy level of the dot as well as of 
the Andreev states. The current-phase relation shows a 
strong modification when an Andreev state intersects the 
Fermi energy and the Josephson junction can be driven 
into the n state by the exchange interaction (Fig. |4][a)). 
A magnetic field applied to the central region induces a 
Zeeman energy of the electrons on the quantum dot and 
a precession of the molecular magnet's magnetisation. 
As a result, an electron on the dot can undergo spin- flip 
scattering and the occupation of the Andreev states is 
redistributed by the magnetic field. The redistribution 
depends on the orientation between the magnetic field 
and magnetisation of the molecular magnet and can lead 
to an increase in the current-phase relation (Fig. p[c)) 
when the electron-like quasiparticles are scattered into 
states below the Fermi energy carrying the current in 
positive direction. The critical current in Fig. W\ shows 
the possibility to observe multiple transitions as a func- 



tion of a gate voltage which changes the energy level of 
the dot. Multiple transitions are also obtained if the 
magnetic field is swept around zero. However, in a more 
realistic system the dependence of the superconducting 
gap on the magnetic field must be taken into account in 
a self-consistent way. 

In principle, the results of the paper can be used in two 
ways. First, the measurement of the physical magnitudes 
of either the energy of Andreev levels, the current-phase 
relation or the critical current, allows to reveal the in- 
ternal structure of the Josephson junction. Different pa- 
rameters are possible to extract from measurements such 
as the strength of the exchange interaction. Second, the 
presence of the molecular magnet offers the possibility to 
manipulate the current through the junction by chang- 
ing the parameters such as the energy level of the dot 
or the magnetic field. This manipulation strongly affects 
the Andreev states w hich a llows possible application in 
quantum computation.'23221 
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